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(ABSTRACT) 


The Thermal Radiation Group, a laboratory in the department of Mechanical Engineering 
at Virginia Polytechnic Institute and State University, is currently working towards the 
development of a new technology for cavity-based radiometers. The radiometer consists of 
a 256-element linear-array thermopile detector mounted on the wall of a mirrored wedge- 
shaped cavity. The objective of this research is to provide analytical and experimental 
characterization of the proposed radiometer. A dynamic end-to-end opto-electrothermal 
model is developed to simulate the performance of the radiometer. Experimental results for 
prototype thermopile detectors are included. Also presented is the concept of the discrete 
Green’s function to characterize the optical scattering of radiant energy in the cavity, along 
with a data-processing algorithm to correct for the scattering. Finally, a parametric study 
of the sensitivity of the discrete Green’s function to uncertainties in the surface properties 
of the cavity is presented. 
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Chapter 1 


Introduction 


Increased industrialization over the last century, while bringing a myriad of benefits to 
society, has greatly impacted the world’s social structure and its ecosystem. Pollution of our 
rivers, oceans, soil, and atmosphere have all increased significantly during this century. 

One area of increasing concern is the ongoing pollution of the atmosphere. Continued emis- 
sions of CFC’s, as well as the burning of fossil fuels, increases the level of CO2 in the 
atmosphere and has led to an increased concern over global warming due to the “greenhouse 
effect.” 

In order to understand the effects of these changes in our climate system, we must first seek 
to understand the climate system itself. For example, it is hypothesized that the greenhouse 
effect causes, through warming of the lower troposphere, increased concentration of water 
vapor in the atmosphere. Because water is itself a greenhouse gas, this raises the possibility 
of “runaway” global heating. However, increased water vapor concentration should increase 
the production of clouds in the atmosphere, which would have an overall cooling effect on 
the climate. 

In an effort to better understand the climate system, the international science community 
has begun a series of programs aimed at increasing our knowledge and understanding of the 


1 



2 


climate through observation of climatic changes and radiation in the atmosphere. 


1.1 Description of the instrument 

The thermopile linear array detector presented here is a new detector concept developed by 
Dr. J.R. Mahan of the Department of Mechanical Engineering at Virginia Tech, and Mr. 
Lawrence W. Langley, president of Vatell Corporation [Mahan and Langley, 1996] , intended 
for use in the Geo- Synchronous Earth Radiation Budget (GERB) instrument. 

The thermopile detector consists of a linear array of 256 pixels, or sensing elements. Each 
pixel is the blackened active junction of a single-junction-pair thermocouple. Each pixel 
is 60 by 60 jum and is separated from the neighboring pixels by a 3-//m laser-etched gap. 
The pixels are assumed to be electrically and thermally independent. The geometry of the 
thermopile detector is given in Figure 1.1 . 

Each element of the linear-array detector consists of a blackened active junction and a 
shielded reference junction. The two junctions are connected electrically by a zinc-antimonide 
bridge and are connected electrically to the external readout circuits by platinum leads in 
each junction. The resulting Seebeck coefficient for the zinc-antimonide/platinum thermo- 
couple is many times greater than that of traditional metal/metal thermocouples. The 
reference junction is in direct thermal contact with an aluminum-nitride substrate, while the 
active junction is separated from the substrate by a layer of parylene, which acts as a thermal 
impedance. A complete description of the thermoelectric performance of the thermopile is 
given by Weckmann [1997]. 

The detector is bonded to one wall of a wedge-shaped cavity, as shown in Figure 1.2. The 
walls of the cavity are highly reflective and have a high specularity ratio, defined as the 
ratio of the specular component of reflectivity to the overall reflectivity. The cavity opens 
to its surroundings through a 60-/un slit directly above the detector. Collimated light enters 
through the slit and strikes the detector. Approximately 90 percent of the radiation is 
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absorbed on direct incidence, with the rest being reflected. However, because of the cavity 
geometry, most of the reflected radiation returns to be eventually absorbed by the detector. 
The portion of the radiation that is reflected specularly is absorbed by the pixel to which it 
was originally incident, while the diffuse portion is generally absorbed by other pixels in the 
array, thereby creating “optical cross-talk.” The elimination of this cross-talk to recover the 
original data is addressed in this thesis. A complete discussion of the optical behavior of the 
cavity and detector is given by Sanchez [1998]. 


1.2 Geo-synchronous Earth Radiation Budget (GERB) 

As stated above, long-term measurements of the earth’s radiant energy budget are needed to 
assist with modeling the climate system and predicting long-term trends. These observations 
need to be satellite based in order to provide sufficient global coverage and sampling. Other 
satellites that carry radiation budget sensors use polar or other low-earth orbits. Using these 
orbits, four satellites can provide coverage of the diurnal cycle with a temporal resolution 
of three hours, which is only marginally adequate. The GERB instrument addresses this 
undersampling by performing measurements from geo-synchronous orbit. This means that 
it always observes the same part of the earth but at a continually changing time of the day. 

The GERB instrument will measure broadband radiances exiting the Earth from the per- 
spective of the ESA’s Meteostat Second Generation (MSG) satellite, with a launch scheduled 
for the year 2000. The instrument will make both shortwave (0.32-4 fin l) and total (0.32-30 
fim) measurements, with the longwave (4-30 /im) component determined from the differ- 
ence between the total and the shortwave. Here “shortwave” refers to the solar spectrum 
and “longwave” to the Earth-emitted spectrum. Each channel has a sampling time of 300 s. 
A more complete description of the GERB instrument is given in Table 1.1 [Mueller, 1997]. 




256-Pixel Linear Array Thermopile Detector 


Figure 1.2: Geometry of the cavity showing placement of the linear-array detector 
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Table 1.1: GERB instrument description [Mueller, 1997] 


Platform 

Meteostat Second Generation 

geo-synchronous (above Odeg, Odeg) 
spin-stabilized (100 rpm) 
launch: 2000 

Spectral Bands 

shortwave 

0.32-4 /m i 


total 

0.32-30 // m 


longwave 

by subtraction 

Short-wave filter 

quartz 


Accuracy (Target) 

shortwave 

1 percent 


longwave 

0.5 percent 

Field of View 

43.7 km 

at sub-satellite point 

Telescope 

3-mirror off-axis optics 



+ 1 fold mirror 

(to minimize polarization) 


+ 1 de-spin mirror 

(to remove the satellite rotation) 

Detector 

256 by 1 

linear array in North-South direction 

Calibration Sources 

thermal infrared/solar 

blackbody /integrating sphere 

Sampling time 

300 s 

(both channels, complete Earth image) 

Mass 

less than 30 kg 

divided into three “boxes” 
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1.3 Green’s functions 

The English mathematician George Green first introduced the idea of potential functions in 
1828 [Green, 1828] . In his essay Green describes a system of bodies S and a material point p 
within the system, and states that the effect, in any given direction, of all the forces acting 
upon the point, may be described by a function that he calls the “potential function arising 
from the system S.” What Green was describing was the response of the system to an impulse. 
In the cited paper, Green also introduces what later became known as Greens functions. 
Green introduced these functions in order to solve the differential equation associated with 
the potential function. Green considered the case where the potential is caused by a single 
charge at an isolated point such that the Laplace equation is satisfied by the potential 
function. By assuming that the potential becomes zero on the boundary of the body and 
increases as 1 jr as the point charge is approached, Green was able to demonstrate how to 
solve for the function [Gray, 1994] . 

Although Green derived his functions to solve problems in electricity and magnetism, the 
underlying principles of his discovery have been extended for use in many areas of science 
and engineering. A particular area that makes extensive use of Green’s functions is that of 
conduction heat transfer, as they are particularly well suited to solve the equation of heat 
diffusion. 

Consider the following three-dimensional nonhomogeneous boundary- value problem of heat 
conduction over region R: 

2 1 . 1 dT(r,t) . . 

V T( r ,t) + -g(r,t) = - — ^ — , *>0, (1.1) 

with boundary conditions on each surface Si defined by 

dT 

h — + hiT = hi T^i = fi (r, t) (1.2) 

and the initial condition 


T(r.t) =F(r), t = 0 . 


(1,3) 



Equations 1.1 through 1.3 can be solved by considering the associated Green’s function 
equation over the same region R: 

V 2 G(r, t\r', t) + j S(r - r')5(t - r) = t > 0 , (1.4) 

k a at 

with the homogeneous boundary condition 

dG 

h — b hi G = 0, t > t . (1.5) 

OTli 

The solution to Equation 1.1 subject to Equations 1.2 and 1.3 is 

T(r,t)= [ G(r,t\r' , t )\ t=0 F (r') dv' 

Jr 

+ |/ dr G(r , t\r', r ) g(r', r ) dv' (1.6) 

+ a f dr^N j G(r,t\r' ,r)| r /= ri ^-/ i (r' ,r) ds[ , 

Jt = o i=1 JSi K i 

where dv' and ds[ are differential volume and surface elements, respectively. 

The physical significance of the Green’s function G(r, t\r' , r) for the three-dimensional prob- 
lem is as follows: It represents the temperature at the location r, at time t, due to an 
instantaneous point source of unit strength, located at the point r', releasing its energy at 
time t = t [Ozisik, 1993]. This definition again illustrates that the Green’s function G 
represents the effect at (r,t) to an impulsive point source at (r', r). 

While the Green’s functions have been used extensively in conduction heat transfer, and to 
some degree in convective heat transfer, their use in radiation heat transfer has remained rel- 
atively rare. Sharma and Cogley [1982] have developed an approach to solving radiative heat 
transfer problems with scattering in a non-gray, plane-parallel environment, using scattering 
functions (Green’s functions). 

With the advent of faster computers, numerical methods of radiation heat transfer analysis, 
such as the Monte-Carlo ray-trace (MCRT) method, make possible a new area of application 
of Green’s functions. Mahan et al. [1998] have proposed the use of a discrete Green’s function 
to represent the scattering in a cavity based on the distribution factor in the MCRT method. 
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1.4 Goals and motivation 

The Thermal Radiation Group, under the guidance of Dr. J.R. Mahan of the Department 
of Mechanical Engineering at Virginia Polytechnic Institute and State University, has been 
involved in the development of optical, radiative, thermal conduction and electronic models 
of spaceborne radiometers for the last twenty-five years. Design and analysis efforts include 
NASA’s ERBE and CERES instruments. 

Over the last three years, the Thermal Radiation Group has been involved in an effort to 
develop a new technology of thermal radiation detectors for use in the next generation of 
spaceborne instruments, including possible use on geo-synchronous earth radiation budget 
(GERB) applications. The technology development has consisted of two major efforts: 

• development of optical and electrothermal numerical models to predict the detector’s 
physical behavior and to aid in the design of the optimal configuration 

• fabrication and testing of actual prototype detectors. 

Towards these ends, a dynamic electrothermal model of the detector was developed by Weck- 
mann [1997] , an optical model of the detector cavity was developed by Sanchez [1998] , and 
an experimental procedure was developed by Barreto [1998] to estimate the electrothermal 
properties of the detector. 

The aims of the current research include: 

• integration of the previous optical and electrothermal models into an end-to-end model 
of the cavity and detector to accurately model and predict the performance of the 
instrument 

• development of a method to recover the incident radiation that has been spread through 
the cavity via reflections, by use of discrete Green’s functions that characterize the 
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scattering of radiant energy due to the cavity effect, as well as performing a parametric 
study of the behavior of the discrete Green’s functions with variations in the physical 
parameters of the cavity 

• testing and characterization of actual physical prototypes of the thermopile radiation 
detectors to determine the actual performance of the detectors. 

A FORTRAN program, Detector. f, is developed to achieve the numerical modeling aims 
listed above. All physical testing of the detectors was performed at the NASA Langley 
Research Center in Hampton, Virginia. 



Chapter 2 


Radiation Detectors 


This chapter attempts to provide the reader with a brief overview of radiation detectors, 
and particularly thermal detectors. Topics include general properties of detectors, including 
responsivity, signal-to-noise ratio (SNR), fundamental noise, etc. The different types of 
thermal detectors are discussed, with particular attention paid to thermopile detectors. It 
should be noted that this chapter is not meant, to be an exhaustive discussion of these topics, 
but rather an introduction to several topics that will arise in later chapters. A large portion 
of this chapter was adapted from material by Lenoble [1993] and Haeffelin [1997] . 


2.1 General properties of radiation detectors 

The input signal, or flux, to a filtered radiation detector is 

-A(A) = «a(A) d> A (A) , (2.1) 

where d>^(A) is the incident monochromatic radiant flux and «a(A) is the monochromatic 
absorptivity of the detector surface at a particular wavelength A. 

The absorptivity of the detector can be enhanced by two means. The first is through the use 
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of a “black” coating with a high absorptivity on the detector surface. The second is through 
a cavity, which causes the incident flux to strike the detector surface multiple times through 
specular (mirror-like) reflections. 

There are two major categories of radiation detectors: thermal detectors and photon detec- 
tors. Thermal detectors transform the absorbed radiant energy into sensible heat, producing 
a temperature difference that can be directly measured, or that can be used to induce an 
electromotive force, or emf. In photon detectors, each photon in the absorbed radiant en- 
ergy is used to excite an electron in the detector material, thereby increasing the energy 
level of the electron. If the energy of the electron is raised sufficiently to cause the electron 
to jump to a higher energy level, the electron is extracted from the detector, producing a 
measurable photoelectric current. 


The responsivity of the detector is defined as the signal output per unit of incident energy, 
and is represented by 


Y Y 
» = — = a — 
$ a X ’ 


where Y is the output yielded by the detector. 

The output of the detector always contains a certain amount of noise. A portion of this noise 
arises from the measurement procedure and can be minimized. The remainder is noise that 
is inherent to the physical properties of the detector, and is known as fundamental noise . 
The output signal of the detector can be expressed as the root-mean-square (rms) value (Y) 
defined 

< 2 - 3 ) 

i = 1 




The signal-to-noise ratio is defined 
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SNR 


(Y) 

({ A y2»l/2 • 


(2.5) 


The noise equivalent power ( NEP ) is the detector input that produces a SNR = 1, and the 
detectivity is given by 


D 


1 

NEP ' 


( 2 , 6 ) 


Detector responsivity is typically wavelength dependent, but different detectors can be char- 
acterized using the normalized detectivity 

D* = D(AAf) 1/2 , (2.7) 


with A being the sensitive area of the detector and A / the measurement bandwidth. 


2.2 Fundamental noise of radiation detectors 

Three types of noise act as fundamental limits on the performance of radiation detectors: 
shot noise, Johnson noise, and temperature noise. 

2.2.1 Shot noise 

Shot noise arises due to the granular structure of the material from which the element 
is manufactured, and affects both the incident radiation and the electrical output of the 
detector. The radiant flux is composed of a beam of n photons striking the detector per unit 
of time, and the electric current density is due to the transport of n electrons per unit time. 
Any measurement procedure of these quantities actually consists of counting n t discrete 
events in a time interval At, the probability of occurrence of which can be represented by a 
Poisson distribution. 
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The average current observed per unit of time is 

' At 


where e is the electron charge of the particle, and the current fluctuation is given by 


(A! 2 ) = <(i, - (i» 2 > . 


(2,9) 


Equation 2.9 can be represented in the frequency domain as 


(Ai 2 ) = 2e(i) A / , 


( 2 , 10 ) 


with A / representing the frequency bandwidth. The voltage fluctuation over a resistor R is 
given by 

(AV 2 ) =2Re{V) Af . (2.11) 


Photons differ from electrons in that the photon energy hu varies with frequency v and the 
radiant energy. For high frequencies (A < 40 fim ) the noise fluctuation is 


(An 2 ) 


( n t) 



1 

exp (hu / kT) — 1 


( 2 , 12 ) 


for photons of frequency v = c/A, where c is the speed of light in the medium, h is the 
Planck constant, and k is the Boltzmann constant. For low frequencies (A > 40 jim) the 
noise fluctuation is 

(2,13) 


(An?) = <»«} 




The corresponding fluctuation in the incident energy is 

(Aw 2 ) = (An 2 ) (hu) 2 = (Ant) hukT , 


(2,14) 


and is known as thermal noise ( J 2 ). At short wavelengths the noise due to radiation is 
dominant, and at long wavelengths it is not. 
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2.2.2 Johnson noise 

Johnson noise is the noise associated with random thermal motion within the resistance of 
the detector. These fluctuations can be regarded as being due to an electromotive force, e/, 
at a frequency /. For a small frequency interval, A/, these fluctuations are given by 

(ej) = EAf . (2.15) 

The constant E has been found to be represented by E = AkT R. where k is the Boltzmann 
constant and R is the electrical resistance of the detector. Equation 2.15 can therefore be 
written 


(ej) = AkTRAf . 


(2.16) 


2.2.3 Temperature noise 


The temperature of a thermal detector fluctuates randomly around a mean value, due to heat 
transfer between the detector and its surroundings. These fluctuations, called temperature 
noise, are expressed by 


(AT 2 ) 


kT 2 
~C t ’ 


(2,17) 


where Ct is the thermal capacitance of the detector. The fluctuation of radiant power is 
given by 

(AXj) = AkT 2 K T Af , (2.18) 


where Kt is the thermal conductance (1 /Rt) of the detector. 
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2.3 Thermal detectors 


The output signal, Y, of a thermal detector is induced by a temperature variation AT caused 
by the absorbed flux. The responsivity of the thermal detector can be expressed as 


K = a 


AT Y 
AT ' 


(2.19) 


The factor AT / X is a physical property that depends only on the thermal characteristics 
of the detector. 


If the temperature of the detector is raised by AT 0 and then allowed to decay without any 
external energy source present, the exponential temperature decay is given by 


AT = AT a exp 



( 2 . 20 ) 


where r = Ct Rt is the thermal time constant of the detector. A smaller time constant 
allows the system to return to equilibrium more rapidly and results in a “quicker” detector 
response. 


If the detector is irradiated with an incident flux X modulated at a frequency / such that 


X =X f exp(2? rj ft) , 


(2.21) 


where j is the imaginary operator, \j — 1, then the temperature difference AT is modulated 
at the same frequency such that 

AT = ATf exp(27r j ft) , (2.22) 


where 


A T f = 


A/ 

K t (1 + 47t 2 f 2 t^) 1 / 2 


(2.23) 


At higher frequencies, ATf falls off as 1//. The maximum operating frequency, / m , is defined 
as 


frt 


1 

2n tt 


(2.24) 



17 


2.3.1 Thermopiles 

Thermopiles operate on the thermoelectric principles on which thermocouples are based. 
The three thermoelectric effects that influence the behavior of the thermocouple are the 
Peltier, Seebeck and Thomson effects [Pollock, 1985]. 

The Seebeck effect, named for German scientist Thomas Johann Seebeck, states that an 
electric current will flow in a closed circuit composed of two dissimilar metals if their junctions 
are maintained at different temperatures. If the circuit is open, an electromotive force, Eab 
is created. This force, termed the relative Seebeck voltage, is directly proportional to the 
temperature difference AT between the two junctions, 

E AB = S AB AT, (2.25) 

where Sab (V/AT) is the Seebeck coefficient. 

The Peltier effect, named for French physicist Jean Charles Peltier, states that when an 
electric current flows across a junction of two dissimilar metals, heat is liberated or absorbed, 
depending on the direction of the current flow relative to the Seebeck voltage associated with 
the junction. The power liberated or absorbed is given by 

P = P AB i, (2.26) 

where i is the electric current and Pab (V) is the relative Peltier coefficient. 

The Thomson effect is due to an electric current flowing through a conductor in which exists 
a temperature gradient. Heat is liberated if the current flow direction coincides with the 
direction of heat flow, otherwise it is absorbed. The power liberated or absorbed is given by 

dT 

P = ^ Tx , ( 2 . 

where dT/dx is the temperature gradient along the conductor and a (T) (V/K) is the Thom- 
son coefficient. 
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The relationship among Equation 2.25, Equation 2.26, and Equation 2.27 was developed by 
Pollock [1985] . Pollock states that “the Seebeck effect is the result of both the Peltier and 
Thomson effects.” The relationship between the three effects is 

= + (2.28) 

For a more thorough treatment of the thermoelectric and thermodynamic properties of ther- 
mopiles, the reader is referred to Weckmann [1997]. 

When used as thermal radiation detectors, one junction of the thermocouple is maintained 
at a reference temperature equal to the ambient temperature, Too, while the other junction 
is typically coated with a black absorber material and allowed to store thermal energy by 
absorbing incident radiation. A thermopile consists of several thermocouples connected in 
series to increase the emf output. For n junction pairs, the output signal is 

Y = AV = nS A B&T. (2.29) 


Huang [1990] presents a formal development of Equation 2.29 as well as excellent methods 
for calibration and signal conversion of thermopiles for temperature measurements. 


The responsivity of a thermopile detector is 


ft 


o: n Sab 

K t { 1 + 4tt 2 /W!) 1 / 2 


(V/W) , 


which, at low frequencies, reduces to 

ft 


a n S A b 
K t 


{V/W) . 


(2.30) 


(2.31) 


2.3.2 Bolometers 

In a bolometer a temperature increase AT due to absorbed heat flux produces a change in 
the resistance R 0 of the detector according to 


A R = R — R a = R a P AT , 


(2.32) 
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Incident Radiation 



Figure 2.1: Schematic diagram of a typical bolometer circuit (adapted from Lenoble [1993]) 


where P (K x ) is the temperature coefficient of the resistance. If the bolometer is connected 
in series with a voltage source, Vj,, and a resistance Rl much larger than R a , as seen in 


Figure 2.1, the voltage V across Rl varies with the temperature according to 

= w Ar . 

Rl 


(2.33) 


The voltage difference A.V is equivalent to the output Y, and thus the responsivity can be 
written as 


V b R 0 P 

R L K t (1 + 4tt 2 / 2 t^) 1 / 2 


(V/W) 


(2.34) 


or at low frequencies 

» = (V/W) . 


(2.35) 


2.3.3 Pyroelectric detectors 

A pyroelectric detector consists of a plate of a crystal that, when heated by absorbed radiant 
flux, produces electric charges on opposite sides of the plate due to a change in average 



20 


position of the electrons in the crystal. This leads to a current in an external circuit if the 
input signal is modulated, as a pyroelectric detector does not respond to an unmodulated 
signal. This is because the initial charge separation drains off through the electrical resistance 
of the crystal material. 


The responsivity of a pyroelectric detector with equivalent electric resistance and capacitance 
of R e and C e , respectively, is given by 


5ft 


a2n f P T A R e 

K t (1 + 47T 2 P Tp) 1 / 2 (1 + 47T 2 P t 2 ) 1 / 2 


(V/W ) , 


(2.36) 


where A is the sensitive detector area, P T is the pyroelectric coefficient, and r e = C e R e is 
the electrical time constant. 

A frequency range, / l/27rr e and / ;§> 1/27T T t exists where the responsivity does not 
depend on the modulation frequency. Pyroelectric detectors are typically operated in this 
frequency range. The responsivity in this frequency range is 

^ ex Pt A R e q Pf A R e 

Kt Tt Ct 


{V/W) . 


(2.37) 


2.4 Thermopile materials 

The following is a brief overview of some of the materials used in thermopile detectors. The 
two areas discussed are materials used in the thermocouple junctions and absorber layers, 
and some of the issues and current technology associated with each. 


2.4.1 Thermocouple junctions 

Two types of thermocouple junction are considered: two thin metal layers, and a thin metal 
layer with a semiconductor material layer. 

An example of a typical thermocouple junction consisting of two thin metal layers is a 
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bismuth-antimony, Bi-Sb, thermocouple. A typical Bi-Sb thermocouple junction has a See- 
beck coefficient of approximately 100 fiV/K [Lenoble, 1993]. 

Recent technological and manufacturing advances have seen an increased use of silicon- 
based thermopile technology. The advantages of using silicon in thermopiles include its 
compatibility with CMOS-standard processes, which results in higher reliability and less 
variation between production runs, and its increased stability with temperature [Schmidt, 
1995]. Another major advantage of silicon-based thermocouples is the ability to alter the 
doping concentration in the material to change the resistivity and thermoelectric power, 
and can be optimized for higher sensitivity /noise ratios. Properly doped silicon can have 
a Seebeck coefficient as high as 1 mV /K , although standard fabrication processes typically 
produce a Seebeck coefficient on the order of 100 fiV jK [Foote et ah, 1998]. 

Recently, new combinations of semiconductor /metal combinations have been under investi- 
gation with the aim of increasing thermopile performance. A thermopile detector comprised 
of zinc-antimony (Zn-Sb) and platinum (Pb) with a Seebeck coefficient of 920 ixV/K [Krieder, 
1994] and a broadband sensitivity of 0.088 [xV/W/m 2 [Mahan et ah, 1998] has been reported. 
Foote et al. [1998] have reported arrays of thermopile detectors comprised of alternating lines 
of Bi-Te and Bi-Sb-Te with reported response times of 99 ms, zero-frequency D* values of 
1.4 x 10 9 cm Hz 1 ' 2 /W and responsivity values of 1100 V/W when viewing a 1000 K black- 
body source. These values represent the best published performance of a thermopile array 
to date. 

2.4.2 Thermopile absorber layers 

The absorber layer coating the hot junction of the thermopile is an essential component 
in the overall performance of the thermopile detector. The absorptivity must be high for 
increased sensitivity, but the mass of the absorber layer must be small for a shorter thermal 
time constant. The deposition of the absorber layer plays a critical role in the development 
of thermopile technology. 
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Three primary methods are available for the deposition of the absorber layer on the ther- 
mopile [Lang et al., 1991] : 

• Thin metallic films. For thick metallic films, the high mobility of electrons in the 
material causes reflection of incident radiation. Conversely, very thin metal films do 
not provide sufficient electrons to interact with the incident light, making the film 
transparent to incident radiation. A critical intermediary thickness exists for which 
the electrons can interact but not move, and thus are able to absorb incident photons. 
A broadband absorption of 0.5 in a very narrow range of film thickness has been 
reported. 

• Highly porous platinum black deposited by a galvanic process. If platinum is deposited 
by electroplating, and a small amount of a heavy metal, such as lead, is added, the 
result is a porous film called “platinum black”. A film of 1.5 mg /cm 2 provides an 
absorptivity of approximately 0.95 at A = 10 gm. 

• Evaporation of a metal in an atmosphere of low-pressure nitrogen. A common pro- 
cedure for adding a porous absorbing layer is the evaporation of gold in a nitrogen 
atmosphere. The resulting layer is called “gold black.” This process is the best of the 
three, providing an absorptivity of approximately 0.95 at A = 10 gm with a film of 250 
gg/cm 2 . 

Some promising “off-the-shelf” coatings include Chemglaze Z-306, ECP-2200, and DeSoto 
Black. When mixed with an additive such as carbon black, these materials exhibit a re- 
flectance of less than 10 percent over a wide spectral range [Smith and Howitt, 1986]. 



Chapter 3 

Testing and Characterization of a 
Thermopile Radiation Detector 

3.1 Design of an output signal amplifier 

As series of prototype detectors was fabricated by the Vatell Corporation and submitted to 
the NASA Langley Research Center for evaluation. Output of these prototype detectors, 
viewing a simulated blackbody at 1100K, was in the range of 50 nV to 50 //V, depending 
on the wavelength band being observed. Available instruments used to read the signal have 
a least count on the order of 1 mV; therefore, the detector output had to be amplified for 
accurate measurements to be taken. 

The design of an appropriate pre-amplifier for this application involved several considera- 
tions. The first of these was determining the gain required to give a measurable signal. A 
gain of 50,000 was decided upon, giving an optical broadband output of approximately 1.5 
V and an optically filtered output of 1-100 mV, depending on the bandpass of the optical 
filters used. The next consideration in the design process was determination of the type of 
amplifier circuit to be used as well as selection of appropriate operational amplifiers (op- 
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Figure 3.1: Schematic circuit diagram of output pre-amplifier 


amps). After several different designs were considered, a multi-stage differential design was 
chosen as being the most appropriate for this application. The first stage of the pre-amplifier 
is a differential stage with a gain of 20, followed by two identical inverting stages with gains 
of 50 each. Once the circuit had been designed, an appropriate op-amp had to be chosen. 
The main consideration in selecting the op- amp is that it be low-noise. This was extremely 
important in this case due to the small amplitude of the detector output signal. Anything 
more than a few millivolts of noise from the amplifier would tend to drown out the detector 
signal and cause a poor signal-to-noise ratio. Several different op-amps were considered and 
tested, and while some had better listed characteristics, the Harris HA-5127 [Harris Semi- 
conductors, 1993] had the best performance when tested under laboratory conditions. A 
schematic diagram of the final circuit is given in Figure 3.1 . 

The frequency response of the pre-amplifier was tested by attaching the input of the pre- 
amplifier to a voltage source of 200 p,V and varying the frequency of the source. The output 
voltage was recorded and the gain calculated for each frequency value. The 3-dB down point 
was found to be approximately 100 Hz. The frequency response characteristics are given in 
Figure 3.2. 


Gain, 



1 10 100 


Log Frequency, Hz 

Figure 3.2: Measured frequency response of output pre-amplifier 




26 


3.1.1 Noise analysis of output pre-amplifier 

A noise analysis of the final amplifier circuit was performed to predict the total level of noise 
that the circuit would contribute to the signal. The sources of noise that were considered 
were voltage and current noise inherent to the op-amps, and Johnson noise, which is thermal 
noise caused by heating of the resistors as current is passed through them. Several other 
possible sources of noise exist, but were considered negligible. Due to the nature of a multi- 
stage amplifier, noise from the first stage of the amplifier will tend to be the dominant noise 
source. This is because feedback in the second and third stages can at least theoretically 
reduce noise in these stages to negligible levels. 

The voltage noise is given by 

N v = e n (l + G ) , (3-1) 

where e n is the voltage noise density of the op-amp and G is the gain of the first stage of 
the pre-amplifier. The current noise of the non-inverting terminal of the op-amp is given by 

N it+ = i n R {+) (l + G) , (3.2) 

where i n is the current noise density of the op-amp and i?( + ) is the electrical resistance seen 
at the non-inverting terminal of the op-amp. Similarly, the current noise of the inverting 
terminal is given by 

Ni~ = « n i?(_)(l + G) , (3-3) 

where R(~) is the electrical resistance seen at the inverting terminal of the op-amp. The 
Johnson noise of each resistor in the first stage is given by 

R n = (4 kTR)* , (3.4) 

where k is the thermal conductivity, T is the temperature, and R is the electrical resistance 
of each resistor. 
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The specified values of the voltage noise density, e n , and the current noise density, i n , for 
the HA-5127 Op-Amp are 3.8 riV/y/ Hz at 100 Hz and 1.7 pA/y/Hz at 100 Hz, respectively 
[Harris Semiconductors, 1993]. 

The calculated noise from the first stage is 3 //V, and the calculated total noise for the 
amplifier is 7 mV. Subsequent testing of the actual amplifier showed a noise contribution 
from the amplifier of approximately 20 mV, roughly three times the predicted value. 


3.2 Testing of the detectors 

The prototype detector provided for testing consisted of five separate detectors aligned in a 
linear array. We were interested in testing for three parameters of these detectors: electrical 
resistance, sensitivity, and spectral response. Specifically, we were not interested in testing 
the dynamic response of these prototypes because they were much larger than the practical 
linear-array detectors. It should be noted that in all ensuing discussions, responsivity refers 
to the output of the detector per unit of radiant power seen by the detector, while sensitivity 
refers to the output of the detector per unit of radiant flux seen by the detector. Of course, 
the sensitivity can be calculated from the responsivity if the area of the detector visible to 
the radiant power is known. 

The electrical resistance of the detectors was an important parameter to determine, both 
because of its effect on the amplifier circuit and to better understand the behavior of the de- 
tectors. This is because a change in electrical resistance with time might imply other changes 
in the properties of the detector. Two tests were utilized to measure the electrical resistance 
of the detectors. The first test involved attaching the positive and negative terminals of each 
detector to a Hewlett-Packard 4285A Precision LCR Meter. Measurements were taken with 
the meter set to a bias of -0.5, 0, 0.5V, and the measurements were averaged over the three 
settings. The second test consisted of using a Hewlett-Packard 4142B Modular DC Source 
Monitor and a METRICS software package. The monitor set a voltage and measured the 
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resultant current across each detector. The software plotted the voltage against current, and 
the slope of the resulting linear fit of the data points gave the electrical resistance of the 
detectors. Results of both tests are given in Table 3.1. 


Table 3.1: Electrical resistance of detectors tested 


Detector 

Electrical resistance, kO 
LCR Meter 

Electrical resistance, kD 
DC Monitor and METRICS 

1 

24.0 

24.12 

2 

9.4 

9.05 

3 

5.4 

5.31 

4 

20.9 

20.92 

5 

20.9 

20.13 


The next parameter of interest was the sensitivity of the detectors. Sensitivity is defined 
as the voltage produced by the detector per unit of flux incident upon the detector. The 
sensitivity of the detectors was tested first with a broadband source, and then with the same 
source viewed through a series of optical filters, to determine if the sensitivity of the detectors 
remained constant over different wavelength bands of incident radiation. 

A benchmark measurement was needed for calibration. This was accomplished using a 
reference detector. A reference detector is typically a similar type of detector with a known 
responsivity or sensitivity. If the two detectors are illuminated with the same amount of 
radiation, then the ratios of the output signals should be the same as the ratios of the 
sensitivities, and from this the sensitivity of the detector of interest can be calculated. 

The detector used as the reference detector was the Oriel IR Pyroelectric Detector, Model 
70124. The responsivity of the detector was 1000 V rms /W and 1920 V v - v jW . 

The test set-up shown in Figure 3.3 was used to test the sensitivity of the detectors. The 
Mikron Blackbody Calibration Source was used as the radiation source for the tests. The 
detector was mounted on an x-y positioning table, positioned such that it lay on the optical 





Figure 3.3: Schematic diagram of testing apparatus 
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axis of the blackbody. The detector surface was mounted perpendicular to and facing the 
source. To ensure that the reference and prototype detectors were positioned in the same 
place with a high degree of repeatability, the focal-plane method was used. A kinematic 
mount was placed in a permanent position in front of and facing the detector with a removable 
microscope objective attached to the mount. The axial position of the detector was adjusted 
until it lay in the focal plane of the objective. This gave repeatable position with 100 
/im precision. The reference scans were taken with a signal chopper in front of the source 
set at 134 Hz, and the output of the Oriel detector running directly into the oscilloscope. 
The prototype detectors were tested with a shutter in front of the source set at 1 Hz, and 
with the output of the detector connected to the input of the pre-amplifier described in 
the previous section, whose output was displayed by an oscilloscope. The prototype and 
reference detectors were tested with no filters, and with three separate filters of different 
passbands. The results of these tests are given in Table 3.2. 

Using a responsivity of the reference detector of 1920 V p - v jW and an area of 19.6 mm 2 , and 
by assuming that the prototype detectors see the same flux as the reference detector, we 
can calculate the sensitivity of the prototype detectors. These results are given in Table 3.3. 
The measured values for broadband sensitivity agree well with the value predicted by the 
end-to-end model of the detector of 870 V/W/cm 2 . 

The spectral response of the sensitivity of the two detectors over a wavelength interval 
2 < A < 12 //rn is shown graphically in Figure 3.4. The sensitivity of each detector decreases 
as the wavelength increases. This is most likely due to decreased absorptivity of the absorber 
coating at longer wavelengths [Smith and Howitt, 1986]. 
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Table 3.2: Signal output for testing of two thermopile detectors 



Mean signal 

Standard deviation 

No filter 



Detector 1 

1.39 V 

0.015 V 

Detector 4 

1.44 V 

0.042 V 

Oriel Detector 

13.85 V 

0.079 V 

Filter K, 2. 0-2. 4 yum 



Detector 1 

86 mV 

4.4 mV 

Detector 4 

86 mV 

6.1 mV 

Oriel Detector 

1.02 V 

0.018 V 

Filter B, 6.548-6.708 yum 



Detector 1 

11 mV 

1.2 mV 

Detector 4 

8.78 mV 

3.7 mV 

Oriel Detector 

197 mV 

0.500 mV 

Filter K, 11.55-12.65 yum 



Detector 1 

2.0 mV 

0.25 mV 

Detector 4 

1.7 mV 

0.13 mV 

Oriel Detector 

114 mV 

2.21 mV 


Table 3.3: Sensitivity values for two prototype thermopile detectors 


Wavelength range, yum 

Detector 1, fiV/W/cm 2 

Detector 4, nV/W/cm? 

Broadband 

852 

882 

2. 0-2. 4 

715 

715 

6.548-6.708 

474 

379 

11.55-12.65 

149 

127 





Chapter 4 


Description of the Analytical Model 


A high-level dynamic opto-electrothermal model was developed in order to accurately sim- 
ulate and predict the behavior of the detector. This model incorporates the Monte-Carlo 
ray-trace (MCRT) method to model the optical behavior of the detector, and a finite-element 
model to characterize the electrothermal behavior of the detector. In addition, an algorithm 
to correct the detector output for optical cross-talk was developed and integrated into the 
model. This end-to-end model was implemented in a FORTRAN program titled Detector. f. 
Figure 4.1 provides a logic flowchart for Detector. f. 


4.1 Optical model 

4.1.1 Detector cavity 

The thermopile detector is mounted in a wedge-shaped cavity, as shown in Figure 1.2. The 
reason for mounting the detector in a cavity is to increase the number of rays absorbed on 
each pixel of the detector, through specular reflections, which leads to enhanced spectral 
flatness. For this reason, the walls of the cavity are highly reflective and highly specular. 
The values for the reflectivity and specularity ratio, defined as the ratio of the specular 
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Figure 4.1: Flowchart for Detector. f 
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component of reflectivity over the overall reflectivity, used in the model are both 0.9. The 
thermopile detector contains an active and a reference junction, with the active junction 
being visible to incident radiation and the reference junction being shielded. The active 
junction requires a high absorptivity and a high specularity ratio. To achieve this, the active 
junction is coated with a black absorber coating such as CHEMGLAZE Z-306. A specularity 
ratio of the detector of less than unity creates optical cross-talk, as some rays are reflected 
diffusely and are absorbed by pixels other than the pixel to which they are initially incident. 
Further discussion of optical cross-talk appears later in the thesis. 

Sanchez [1998] performed a study of the optical characteristics of the cavity and determined 
that a wedge angle of 45 deg provides the most uniform distribution of rays across the surface 
of the detector and the least optical cross-talk. For these reasons this was the angle used in 
the current end-to-end model. 

4.1.2 Monte-Carlo ray-trace model 

The optical behavior of the cavity was modeled using the Monte-Carlo ray-trace (MCRT) 
method. A FORTRAN program titled GERB, developed by Sanchez [1998], was modified 
and used as a subroutine in Detector. f. The program reads an input file that can be modified 
by the user, and that specifies the radiometric flux initially incident to each of the 256 pixels 
in the linear array. It then uses the MCRT method to calculate a modified distribution factor 
(Dh). The original distribution factor defined (Dy ) is defined as the number of rays emitted 
from surface i that are absorbed by surface j over the total number of rays emitted from 
surface i, for each pixel [Mahan and Eskin, 1984], The modified distribution factor (ZA) is 
defined as the number of rays incident to surface j that are absorbed by surface i over the 
total number of rays incident to surface j. These values are then used to calculate the total 
flux absorbed by each pixel. 
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4.1.3 Optical cross-talk 

As stated above, the diffuse component of reflectivity of the surface of the active junction 
and the cavity walls causes energy to be absorbed by pixels on the detector array other than 
the pixel on which it was initially incident. This is termed optical cross-talk. It is possible 
to eliminate this cross-talk through post-processing and to recover the values of the original 
incident energy. The post-processing algorithm is based on the “discrete” Green’s function. 

Mahan et al. [1998] have shown that the energy absorbed, X*, on pixels i = 1 . . . n can be 
represented as the sum of the absorbed portion of the incident energy, X“ bs , and the energy 
absorbed due to reflections, X” e ^, 

Xi = Xf s + X[ e/ . (4.1) 

The energy absorbed on initial incidence can be expressed as the product of the absorptivity 
of pixel i, a*, and the energy originally incident to i, The energy absorbed due to 
reflections can be expressed as the product of the discrete Green’s function Gjj and the 
radiation initially incident to pixel j, d>j, for pixels j = 1 . . . n. The discrete Green’s function 

is defined as the radiation initially incident to surface j that is reflected and ultimately 

absorbed by surface i. Equation 4.1 can then be rewritten 

Xi = + GijQj . (4.2) 

Introducing the Kronecker delta function, defined 

■ (43) 

Equation 4.2 can rewritten 

^ = {aAj + Gij)$j . (4.4) 


We can utilize the properties of the MCRT method to deconvolve the output. The distribu- 



(4,5) 


tion factor D' t] in the current application is estimated 

, Na 

ry *J 

ij ~ N, ’ 

where Nj is the number of rays initially incident to surface j and N tJ is the number of rays 
initially incident to surface j that are absorbed by surface i. The matrix D\ - is obtained by 
illuminating each pixel j separately with a fixed number of rays (typically on the order of 
10 6 ) and recording the number of rays absorbed in each pixel i. The result is a n-by-n matrix 
of D[a values, where n is the number of pixels in the linear array, that completely describes 
the distribution of rays between the aperture and the linear array detector (it should be 
noted that this matrix does not represent rays absorbed by the cavity walls). Using the £>h 
matrix, the relationship between the energy originally incident to the pixels (dq) and the 
energy absorbed by the pixels (Xj) is simply 

[X f ] = (£>',][#,] . (4,6) 

We see that Equation 4.6 is in the same form as Equation 4.4 and that D\ - is identical to 
the quantity arfij + Gij. 

Since the matrix is only dependent on the surface properties and geometry of the cavity, 
we can recover the original input to the cavity as long as we can determine the energy 
absorbed on each pixel. This is accomplished by a simple matrix inversion, that is 

w = py-w , (4-r) 

where d>' is the estimated flux initially incident to pixel j. 

The value for the absorbed energy on each pixel, Xj, is related to the pixel output voltage 
according to 
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where 9ft is the responsivity of the detector in V/W . Some electrical noise will be added to 
the signal before the absorbed energy can be determined. This will compromise the accuracy 
of the recovered data, and is the subject of a study later in the thesis. 


4.2 Thermoelectric model 

In order to predict the output of the detector, it is necessary to compute the temperature 
distribution in each pixel, which consists of a thermopile thermal radiation detector. This is 
achieved by creating a thermoelectric model of the detector using a commercial finite element 
software package, ALGOR [ALGOR, 1993]. 

The proposed sensor design is a linear-array thermal radiation thermopile detector. The sen- 
sor uses a differential thermopile arrangement to detect incident radiant heat flux (energy 
per unit time per unit area). The detector is composed of a pattern of thermocouple junction 
pairs, each containing an active junction and a reference junction, electrically connected in 
series. Referring to Figure 4.2, the active junction consists of an absorber layer, a layer of 
zinc-antimonide, a layer of platinum, and a layer of parylene bonded to an aluminum-nitride 
substrate. The reference junction consists of a layer of zinc-antimonide and a layer of plat- 
inum bonded to the same aluminum-nitride substrate. The two junctions are electrically 
connected by a zinc-antimonide bridge, that connects the zinc-antimonide layer in the two 
junctions. As heat flux is absorbed by the two junctions, the reference junction is able to 
dissipate heat directly to the substrate and thus more rapidly than the active junction, due 
to the insulating layer under the active junction. This causes the temperature of the active 
junction to rise more rapidly than that of the reference junction, and so a temperature dif- 
ference is produced. The two junctions produce voltages proportional to their temperatures 
but with opposing polarities. The net result is a positive voltage, equaling the algebraic sum 
of the two, that is proportional to the temperature difference and, at steady state, propor- 
tional to the radiant heat flux. The output voltage is related to the temperature difference 



39 


25.4 pm 



Absorber Layer, 10 jim 
Zinc-Antimonide, 1 pm 
Platinum, 1 pm 
Parylene, 25.4 pm 
Aluminum-Nitride substrate 


Figure 4.2: Two-dimensional geometry of a single thermocouple pair 


between the two junctions by 

V = S ab AT , (4.9) 

where Sab is the Seebeck coefficient having a nominal value of 920 gV/K for the platinum/zinc- 
antimony combination used in the detector. 

The two-dimensional geometry of each pixel, or thermocouple pair, is represented in Fig- 
ure 4.2. The nominal material properties are given in Table 4.1. For more information on 
the thermoelectric characteristics of the detector, the reader is referred to Weckmann [1997] . 

Table 4.1: Nominal material properties 



Mass Density 
(kg/m 3 ) 

Specific Heat 
(J/kgK) 

Conductivity 

(W/mK) 

Platinum 

21450 

133 

71.6 

Zinc Antimony 

6880 

200 

60 

Aluminum Nitride 

3260 

800 

165 

Parylene 

1289 

712 

0.084 

Chemglaze Z-306 

1400 

669 

0.209 
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4.2.1 General assumptions and approximations 

The following assumptions and approximations were made in the thermoelectric model: 

• Conduction and radiation are the only modes of heat transfer present in the analysis, 
as the detector and cavity operate in a vacuum environment. 

• Heat transfer via radiation between the detector and cavity walls is considered negligi- 
ble since the walls and substrate are maintained at a constant temperature of 311 K, 
which is within a fraction of a milhkelvin of the temperature throughout the pixel. 

• No significant variation in the material properties with temperature occurs, as the 
temperature differences and variations are small. 

• No electrical current passes through the thermocouple junctions because the signal- 
conditioning electronics used to monitor the voltage output has an effectively infinite 
impedance. 

• Each thermocouple junction is isolated from the neighboring thermocouple junctions 
by a 3-/jrn vacuum gap. A heat flux boundary condition is applied to the active 
junction, which represents the net absorbed portion of the incident radiation from the 
scene being viewed. The substrate is maintained at 311 K, and represents a constant 
temperature boundary condition at the bottom of the thermocouple junction. All other 
surfaces of the junction pair are insulated, as indicated in Figure 4.3. 


4.2.2 Finite element model 

A two-dimensional finite element model was created to determine the temperature distribu- 
tion in the thermocouple junction. A mesh consisting of 4007 nodes was used, with a constant 
heat flux applied to the absorber layer and with the bottom of the junction maintained at 
a constant uniform temperature. A 10-|Um thick portion of the substrate was included in 
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Heat Flux 1 W/m 2 

UUUUUUi 


Temperature = 31 1 K 

Figure 4.3: Boundary conditions for two-dimensional heat transfer analysis of a thermocouple 
junction 

the model, to test for possible conduction through the substrate between the active and 
reference junctions. A variable shift 6 = T — 311 was applied because of the extremely 
small temperature gradients in the junction (on the order of 10 -4 K). The software used, 
ALGOR, supplies a maximum of four significant figures, so an applied temperature of 311 K 
would cause the actual temperature variation to be submerged in round-off error. Using the 
variable shift, the base temperature becomes 0 K and the temperature range in the junction 
is then on the order of 10~ 4 K. 


4.3 Noise analysis 

In order to accurately model a real instrument, it is necessary to simulate the Gaussian noise 
associated with the detector and the electronic components of the instrument. This noise 
may be thought of as a pertubation to the signal. 

To simulate the noise, we begin with the assumption that the noise associated with the signal 
has a Gaussian, or normal, probability distribution function (pdf) given by 

-(x-vr 

2( cr ) 2 . 



pdf ( x ) = 


<7\/27T 


exp 
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The cumulative distribution function ( cdf ) is defined as 


~df{z) = [ pdf (x) 

J — OO 


(4.11) 


and gives the probability that a normally distributed random variable Z will be less than or 
equal to z. Because noise has a positive value between 0 and oo , the lower bound of the 
integral can be changed to 0 and Equation 4.11 can be written 


cdf (z) = 


rz 1 \-(x-/j,) 2 ] 

^^ ex P v, dx . 


I o ay/2n 


(4.12) 


Evaluation of this integral gives 


cdm = -\ e rf F-erf ^ 


(4.13) 


The error function can be well approximated by 


erf{x) = y? S'- 1 *” ui 


n + 1) 


(4.14) 


Combining Equations 4.13 and 4.14 gives 


A 2n+1 


1 °° 

Cdf(Z) = “V? «! (2n + 1) + X § (_1) " n\ (2n 4- 1) ’ 


j^2n+l 


(4.15) 


where 


A _ V2 (-z + n) 

2 a 


The cumulative distribution function has a property such that “if X is a random number, 
and we have the cdf, F x (X), and if F x (X) is continuous, the random variable Y produced by 
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the transformation Y = F x (X) has a uniform distribution over the interval (0,1)” [Gibbons 
and Chakraborti, 1992], This property allows us to generate an equally distributed random 
number between 0 and 1, assign this value to cdf (x), and solve for x. This generates 
a normally distributed random variable x from a uniformly distributed random number 
generator. 

This principle was applied to the end-to-end model to add noise to the output signal of the 
detector. A value for the standard deviation, <r, is specified in the program, and the value for 
the output signal of the thermocouple junction pair is assigned to the mean, /x. A random 
number between zero and unity is generated, and Equation 4.15 is used to calculate a new 
value for the signal, which lies on a normal distribution around the mean, or original signal 
output. This procedure is repeated for each thermocouple junction pair. 

To test the validity of this procedure, a value of unity was assigned to /x and a value of 
0.1 was assigned to a and n = 10,000 noise values were calculated using the subroutine 
noise2.f. A histogram of these values is shown in Figure 4.4. The normal distribution curve 
is well represented by the histogram, with some small error in the wings due to truncating 
the infinite series expansion of Equation 4.15. 



Density 


Histogram for output of noise2.f 

(p= 1 . 0,0 = 0.1, n = 10000) 



0.7 0.8 0.9 1.0 1.1 1.2 1 

Perturbed signal 


Figure 4.4: Histogram of noise distribution 



Chapter 5 

Results and Discussion 

5.1 Finite element analysis 

5.1.1 Thermal cross-talk 

One of the assumptions made in formulating the thermoelectric model of the detector is 
that each of the 256 individual thermocouple junctions (pixels) is thermally isolated from 
its neighbors by a 3-//m gap that penetrates approximately 10 [mi into the aluminum- 
nitride substrate. To test this assumption, a thermal analysis of two neighboring pixels was 
performed. It was assumed that radiative exchange between the two pixels was negligible, 
and the only mode of thermal interaction was conduction through the detector’s substrate. 
Furthermore, it was assumed that all interactions would most likely occur between the active 
junctions, as the reference junctions are maintained at the temperature of the substrate. A 
second two-dimensional model was created, this time in the y-z plane of the thermocouple- 
pair, at x = 0, with 100 //m of the substrate included to allow for conduction through it. 
Symmetry was assumed to justify an insulated boundary condition to the two sides of the 
model. The bottom of the included portion of the substrate was given a uniform temperature 
boundary condition of 0 K. One pixel was given a flux boundary condition of 1 W/m 2 while 
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Figure 5.1: Two-dimensional temperature distribution of two neighboring active thermocou- 
ple junctions in the y-z plane at x = 0 

the other was allowed to thermally interact with the cavity at 6 = — 313/i = OK through 

thermal radiation. The steady-state temperature distribution of the two pixels is shown in 
Figure 5.1. From the figure it is clear that conductive heat transfer between the pixels is 
negligible. This is most likely due to the thermal impedance layer under each pixel, which 
acts not only to create a temperature difference between the active and reference junctions 
of each thermocouple-pair, but also to isolate them from their neighbors. 
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5.1.2 Two-dimensional model 

In order to determine the output of the detector, it is necessary to accurately determine 
the temperature distribution in the active and reference junctions of the thermocouple junc- 
tion pair. The question arises of which temperature to use in Equation 4.9, however, as 
the temperature is not uniform throughout the active junction. Weckmann [1997] argues 
convincingly that the output of the thermocouple junction should be dependent only on 
the highest temperature in the active junction due to the law of intermediate temperatures, 
and that the difference between this highest temperature and that of the reference junction 
should be used with the Seebeck coefficient to determine the output. This is the approach 
used here. 

In order to determine where in the active junction the highest temperature occurs, the two- 
dimensional finite element model of the thermocouple junction was used under steady-state 
conditions with the active junction illuminated with a heat flux of 1 W/m 2 . The resulting 
steady-state temperature distribution for the thermocouple junction is shown in Figure 5.2. 

It was found from the temperature distribution that the temperature in the active junction 
does not vary in the y-direction, but does vary somewhat in the x-direction. Figure 5.3 shows 
the temperature in the active junction from x = 1 to x = 60 \xm. The highest temperature 
in the active junction with uniform heating is at x = 0 /rra, and the temperature at this 
point is used in all subsequent cases of uniform illumination to determine the voltage output 
of the thermocouple junction. 

Another property of the thermocouple junction that is of interest is the spatial sensitivity 
of the detector to a “limited” source. In other words, the active junction is illuminated by 
a limited source at different locations to determine if the detector response is sensitive to 
position. The two-dimensional finite element model was once again utilized to study this 
question. Due to the symmetry in the geometry of the detector, there should be no variation 
in the response of the detector with distance from gap to gap (z-direction), and so only 
variations in the distance from edge to edge parallel with the gap between the active and 
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Figure 5.2: Two-dimensional steady-state temperature distribution for a thermocouple junc- 
tion with the active junction illuminated by a heat flux of 1 W/m 2 
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Figure 5.3: Active junction maximum temperature as a function of horizontal distance from 
edge 
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x, (xm 

Figure 5.4: Highest temperature of active junction at it is illuminated at location x by a 
5 -pm wide l-W/m 2 source 

reference junctions (x-direction) were considered. A 5 -fim wide heat flux source of l-W/m 2 
strength was applied at 12 locations, each 5 / im apart, across the 60-/irn wide active junction. 
The resulting highest temperature in the active junction for each of the 12 cases is plotted 
in Figure 5.4. The highest temperature in the active junction appears to increase when the 
active junction is illuminated at a position closer to the bridge that separates the active 
and reference junctions. A much stronger effect but having the same trend was observed 
experimentally using laser heating [Kist, 1998]. This behavior could be used to aid in future 
designs of the thermopile geometry. 
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5.1.3 Determination of the thermocouple junction transfer func- 
tion 

The two-dimensional finite element model permits calculation of the temperature distribution 
in the thermocouple junction and the associated voltage output. However, due to the nature 
of the finite element package used, this is not a practical means of determining the output 
of the thermocouple junction within the larger program, Detector. f. Therefore, a transfer 
function for the thermocouple that could be inserted into Detector. f is needed. At this 
stage, for simplicity, only the zero-frequency transfer function of the thermocouple junction 
was considered. With this simplification, the transfer function is a function of the absorbed 
flux alone and can be determined without significant difficulty. 

The two-dimensional finite element model was run for ten cases, for absorbed flux values 
X = 0.5, 0.6, . . . , l.bW/rri 2 . The resulting highest temperature difference between the active 
and reference junctions is shown as a function of absorbed flux in Figure 5.5. The resulting 
plot is linear, and applying a linear fit over the values results in the desired zero-frequency 
transfer function, given by 

AT = 9 x 10“ 12 + 10“ 5 X . (5.1) 

5.2 Optical cross-talk 

The proposed detector design makes use of the optical cavity to increase the relative absorp- 
tivity of the linear-array of thermocouple pairs. One unfortunate side-effect of the cavity, 
however, is that a small fraction of the incident radiation is reflected and absorbed by pixels 
other than those of initial incidence. This can lead to a small error in the output of each 
pixel. 

To determine what effect optical cross-talk has on the instrument, the end-to-end model, 
Detector. f, was used. An input flux of 1 W/m 2 was specified on pixels 101 through 156, 



Temperature Difference, pK 
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Figure 5.5: Temperature difference between the active and reference junctions of a thermo- 
couple junction as a function of absorbed flux 
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with the remaining pixels not being illuminated. The output of each pixel was recorded 
and converted to a flux value, using Equation 5.1. The optical cross-talk on each of these 
pixels was calculated by subtracting the amount of energy initially incident to each pixel (in 
this case, 0.9 W/m 2 , corresponding to an absorptivity of 0.9 for the illuminated pixels). The 
optical cross-talk, prior to data post-processing, is shown in Figure 5.6. The greatest fraction 
of the reflected energy remains in pixels 101, . . . ,156, contributing to the enhanced relative 
absorptivity. However, some spread of energy occurs to pixels not originally illuminated, as 
can be seen in the rounding of the edges of the original “square” input, and in the small 
wings that have developed outside the original illuminated region. 

The discrete Green’s function was utilized to recover the original input energy distribution. 
Recall from Equation 4.7 that the original energy incident to each pixel can be recovered 
by inverting the discrete Green’s function matrix and multiplying it by the output for each 
pixel. This method was employed to restore the scattered radiation to the pixel upon which 
it was initially incident. Figure 5.6 shows the distribution of scattered radiation over pixels 
i = 1, . . . , 256. The recovered energy distribution over pixels i = 1, . . . , 256 is shown in 
Figure 5.7. The scattered energy has been put back into the originally illuminated pixels 
and the original shape of the incident distribution across the array has been recovered. An 
important result from Figure 5.7 is that the algorithm has corrected the output 
for the radiation absorbed by the cavity walls and lost out the aperture! The 
maximum value of A, in Figure 5.6 is slightly higher than 0.05, which indicates approximately 
4 to 5 percent of the radiation has been lost to the walls or the aperture. However, the value 
of <f>j' in Figure 5.7 for the illuminated pixels is 0.1, which indicates that all the radiant 
energy has been recovered. 
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5.3 Study of detector noise 

5.3.1 Effect of data deconvolution on signal-to-noise ratio 

It has been shown that the original incident energy can be recovered almost perfectly from the 
absorbed energy distribution once the discrete Green’s function is known. A real instrument, 
however, would have to recover the original data from a signal with some amount of noise 
added by the detector and associated electronics. One concern is how the inversion technique 
using the discrete Green’s function, which is essentially a differentiation process, amplifies 
the noise level of the signal. 

To investigate the possible effect of the use of the discrete Green’s function inversion tech- 
nique on a signal with Gaussian noise, the program Detector. f was modified to add a noise 
component with a Gaussian distribution to the output signal. The program was run for 
three different cases, each representing a different signal-to-noise ratio. The results of these 
three cases are given in Table 5.1. 

For the values in Table 5.1, the mean signal was calculated by taking the mean of the signal 
values (Y) for pixels i = 101, . . . , 156 and noise was calculated by considering the variance 
(Y 2 ) of the values over this interval. The same method was applied to the recovered flux <&' . 

For all three cases presented, the signal-to-noise ratio decreases slightly after the inversion 
technique is applied. The decrease is approximately 1 percent for all three cases. 

Figure 5.8 is a plot of the output signal of the detector with Gaussian noise added, before 
the discrete Green’s function inversion is applied, for the case of a signal-to-noise ratio of 
approximately 10. Figure 5.9 is a plot of the same signal after the inversion technique is 
applied. 
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Figure 5.8: Thermocouple junction output signal with Gaussian noise added, signal-to-noise 
ratio of 10 
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Table 5.1: Signal and noise information for output signal of the thermocouple detector prior 
to and following data processing to recover original signal 



Mean Signal 

Noise 

Signal-to-noise 

% difference 

Case 1 

output signal Y, nV 

8.7454 

3.8211 

2.289 


recovered flux , W / m 2 

0.99572 

0.43934 

2.266 

+0.9749 

Case 2 

output signal Y , nV 

8.7442 

1.9097 

4.579 


recovered flux +' , W/m 2 

0.99558 

0.21959 

4.534 

+0.9832 

Case 3 

output signal Y, nV 

8.7611 

.95405 

9.183 


recovered flux +' , W/m 2 

0.99750 

0.02618 

9.091 

+0.9990 


% difference = 


SNRy - SNR&. 
(SNR) 


5.4 Parametric study based on the discrete Green’s 
function 

As stated in Chapter 4, the discrete Green’s function, Gij, is the fraction of radiant energy 
initially incident to pixel j that is ultimately absorbed by pixel i. This is of interest since it 
provides a means to measure the scattering in a cavity with the MCRT method. Previously 
only the total absorbed flux was quantified with this method, through use of the modified 
distribution factor. Although the discrete Green’s function is closely related to the modified 
distribution factor, there is a subtle difference, and it is of some value to separate the two. 

Of some interest is the relationship between the discrete Green’s function for a particular 
cavity and the physical parameters that govern the scattering in that cavity. In the MCRT 
environment, scattering is quantified by two parameters: the absorptivity, a, defined as the 
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fraction of energy incident to a surface that is absorbed; and by the specularity ratio, r s , 
defined as the fraction of energy reflected from a surface that is reflected specularly. It would 
be of some value to study the effect of these parameters on the discrete Green’s function, with 
the goal of determining a function to relate them. This would provide a method to determine 
the distribution of energy in a cavity for any radiant input if the physical parameters are 
known. 

A parametric study was performed to study the sensitivity of the discrete Green’s functions 
to variations in the absorptivity and specularity ratio of the detector surface. The wedge- 
shaped cavity used in the end-to-end detector model considered elsewhere in this thesis 
was used in this study also. The program Detector. f was modified such that only one 
pixel, j = 128, was illuminated, and the absorbed energy on each pixel i = 1, . . . , 256 was 
recorded and written to a data file. The program was run for twenty-five cases, corresponding 
to absorptivity values a = 0.5, 0.6, . . . , 0.9 and specularity ratios r s = 0.5, 0.6, . . . , 0.9. 

The discrete Green’s function, Gi % \ 2 8, for & = 0.5 and specularity ratios r s = 0.5, . . . , 0.9 
is shown in Figure 5.10. Only the Green’s functions for i = 100, . . . , 156 were plotted, as 
all values outside this range were essentially zero. In the plot, it is somewhat difficult to 
ascertain the variation of Gi t 128 over i as the values at i = 128 tend to dwarf the values of 
the neighboring pixels. To remedy this, Figure 5.11 presents the same results, but with the 
vertical scale adjusted to display the values of the pixels neighboring pixel i = 128. The 
same vertical and horizontal scale is used for Figure 5.12 through Figure 5.15, which provide 
plots of Gi t 128 for a = 0.6, . . . , 0.9 and the same range of specularity ratios as in Figure 5.10. 

Analysis of Figure 5.11 through 5.15 yields two important facts. First, the values of Gi.us 
decrease as the specularity ratio and the absorptivity increase, and the distribution on i 
becomes more tightly centered around i = 128. This is of course due to there being less energy 
scattered as the absorptivity increases and less being scattered diffusely as the specularity 
ratio increases. More interestingly, the shape of the distribution over i remains similar and 
seems to be scaled with the specularity ratio values, at least over the range of absorptivity 
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values tested. This suggests that there may exist a relatively simple function that relates 
these parameters that could be obtained from a multi-variable regression. Surface plots of 
G\ 27,128 and (J 120, 128 f° r & = 0.5,... ,0.9 and r s = 0.5,... ,0.9 are given in Figures 5.16 
and 5.17, respectively, to illustrate the relationship between the discrete Green’s function, 
absorptivity, and specularity ratio for two individual pixels. The discrete Green’s function is 
shown to be decreasing as both the specularity ratio and the absorptivity increase. This is 
true for both pixel i = 127 (Figure 5.16) and pixel i = 120 (Figure 5.17), with the maximum 
value for the discrete Green’s function being smaller for i = 120, since the pixel is further 
from the source of the impulse. Due to the symmetrical nature of the scattered radiation 
distribution, it is assumed that pixels % = 129, . . . , 256 would have discrete Green’s functions 
similar to pixels i = 1, . . . , 127 

Of particular interest is the discrete Green’s function where i = j, that is the percentage 
of reflected energy initially incident to a pixel that is absorbed by that same pixel. This 
represents the true cavity effect, which is an increase in apparent absorptivity of a given pixel 
due to specular reflections. A surface plot of Gi 28 ,i 28 for a = 0.5, . . . , 0.9 and r s = 0.5, . . . , 0.9 
is shown in Figure 5.18. It is interesting to note that, unlike the discrete Green’s functions 
for i j, which decrease as both the absorptivity and specularity ratio increase, the Green’s 
function for i — j decreases as absorptivity increases but also increases with the specularity 
ratio, which is opposite the cases where i ^ j. This can be explained from the fact that, 
as the specularity ratio increases, less radiation is reflected diffusely, causing more of the 
reflected radiation to be absorbed by the pixel to which it was initially incident, and less to 
be absorbed by the neighboring pixels. 

A logical extension of this parametric study is the development of a function that relates 
the discrete Green’s function to the absorptivity and specularity ratio of the detector surface 
for all the pixels in the linear-array detector. This can be done in the MCRT environment 
by illuminating individual pixels and recording the distribution factors as defined in this 
thesis for each pixel, over a range of values for the specularity ratio and absorptivity. A 
multi- variable regression could then be performed to determine G t j for i = 1, . . . , 256 and 
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Figure 5.10: Discrete Green’s functions, G'^ 128 , for a = 0.5, r s = 0.5, ... ,0.9 

j = 1, , 256. Once a relational function is determined, it could be used to characterize the 
absorptivity and specularity ratio of an unknown detector surface. By actually illuminating 
an individual pixel, using a laser beam, and measuring the absorbed energy in each pixel, 
the discrete Green’s function for each pixel could be estimated, and the surface properties 
backed out of the functional relationship. 
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Figure 5.16: Surface plot of the discrete Green’s function, G\ 27 , 128 ; of pixel i = 127 for an 
impulse on pixel j = 128, for absorptivity and specularity ratios 0.5, . . . , 0.9 
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Figure 5.17: Surface plot of the discrete Green’s function, G 120 , 128 ) of pixel i = 120 for an 
impulse on pixel j = 128, for absorptivity and specularity ratios 0.5, . . . , 0.9 
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Figure 5.18: Surface plot of the discrete Green’s function, Gi 28 ,i 28 , of pixel i = 128 for an 
impulse on pixel j = 128, for absorptivity and specularity ratios 0.5, . . . , 0.9 


Chapter 6 

Conclusions and Recommendations 

6.0.1 Conclusions 

The following conclusions can be drawn from the research presented in this thesis: 

• the two prototype detectors that were tested in the course of this effort had sensitivity 
values of 883 and 826 ^V/W/cm 2 , respectively 

• the end-to-end model of the instrument predicted a sensitivity of 870 /iV/W/ cm 2 for 
a 1 W/m 2 incident flux, which agreed well with experimental results 

• a zero-frequency transfer function was found for the thermopile detector, which pre- 
dicted a linear response when the detector was in equilibrium with a steady radiant 
flux input 

• No thermal cross-talk is predicted between adjacent pixels in the linear array detector 
due to conduction through the detector substrate 

• optical cross-talk was compensated for through use of the discrete Green’s functions 
for an input without noise, and a signal with Gaussian noise was recovered with a one 
percent decrease in the signal-to-noise ratio 
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• a parametric study of the discrete Green’s function, Gij, was performed, and it was 
determined that decreases as both the specularity ratio and absorptivity increase 
for the case of i ^ j, and G l3 decreases as the absorptivity increases but increases as 
the specularity ratio increases for the case i = j. 

6.0.2 Recommendations for further research 

Suggestions for areas of further research include: 

• further studies should be made into the nature and behavior of the discrete Green’s 
functions in radiometric cavities containing detector arrays for various geometries and 
surface properties 

• different geometries and materials should be developed to increase the performance of 
the thermopile detector, including but not limited to absorber materials 

• prototype thermopile detectors should be further tested to provide a more complete 
characterization of their performance 

• diffraction effects should be integrated into the end-to-end model of the cavity radiome- 
ter 

• the optical cavity model should be used in conjunction with the discrete Green’s func- 
tion to measure the values of absorptivity and specularity ratio. 



References 


ALGOR, 1993, Standard User’s Manual. 

Barreto, J., 1998, “Experimental Design for Estimating Electro-Thermophysical Properties 
of a Thermopile Thermal Radiation Detector,” M.S. Thesis, Virginia Polytechnic Institute 
and State University. 

Foote, M. C., Jones, E. W., and Caillat, T., 1998, “Uncooled Thermopile Infrared Detector 
Linear Arrays with Detectivity Greater than 10 9 cm Hz 1 / 2 /W IEEE Transactions on 
Electron Devices , Vol. 45, No. 9, pp. 1896-1901. 

Gibbons, J. and Chakraborti, S., 1992, Nonparametric Statistical Inference, Third Edition, 
Mercel Dekker Inc., New York, NY. 

Gray, J., 1994, “Green and Green’s Functions,” The Mathematical Intelligencer, Vol. 16, 
No. 1, pp. 45-47. 

Green, G., 1828, “An Essay on the Application of Mathematical Analysis to the Theories of 
Mathematics and Magnetism,” Nottingham, England. 

Haeffelin, M., 1997, “Introduction to Atmospheric Science,” ME 4984, Virginia Polytechnic 
Institute and State University. 

Harris Semiconductors, 1993, “HA-5127 Ulta-Low Noise Precision Operational Amplifier,” 
Harris Corporation. 


72 



73 


Huang, B., 1990, “A Precise Measurement of Temperature Difference Using Thermopiles,” 
Experimental Thermal and Fluid Science , Vol. 3, No. 3, pp. 265-271. 

Kist, E., 1998, Personal communication. 

Krieder, K., 1994, “Intermetallic Thermocouples,” United States Patent, Patent Number 
5,356,485. 

Lang, W., Kuhl, K., and Sandmaier, H., 1991, “Absorbing Layers for Thermal Infrared 
Detectors,” TRANSDUCERS ’91. 1991 International Conference on Solid-State Sensors 
and Actuators, pp. 631-4. 

Lenoble, J., 1993, Atmospheric Radiative Transfer , A. DEEPACK Publishing, Hampton, 
Virginia. 

Mahan, J. R. and Eskin, L. D., 1984, “The Radiation Distribution Factor - Its Calculation 
Using Monte Carlo and Example of Its Application,” First UK National Heat Transfer 
Conference, pp. 1001-1012, July 4-6, Leeds, Yorkshire, England. 

Mahan, J. R. and Langley, L. W., 1996, “The Geo-Synchronous Earth Radiation Budget In- 
strument: A Thermopile Linear- Array Thermal Radiation Detector,” Proposal submitted 
to NASA, Hampton, VA. 

Mahan, J. R., Weckmann, S., Sanchez, M. C., Sorensen, I. J., Coffey, K. L., Kist, E. H., and 
Nelson, E. L., 1998, “Optical and Electrothermal Design of a Linear- Array Thermopile 
Detector for Geostationary Earth Radiation Budget Applications,” SPIE Proceedings , In- 
ternational Symposium on Remote Sensing, Vol. 3498. 

Mueller, J., 1997, “Geostationary Earth Radiation Budget (GERB) Instrument Calibration 
Plans,” Advances in Space Research , Vol. 19, No. 9, pp. 1307-16. 

Ozisik, M. N., 1993, Heat Conduction, Second Edition, John Wiley & Sons, New York. 



74 


Pollock, D. D., 1985, Thermoelectricity, Theory, Thermometry, Tool , ASTM Special Tech- 
nical Publication, American Society for Testing and Materials, Philadelphia, PA. 

Sanchez, M. C., 1998, “Optical Analysis of a Linear- Array Thermal Radiation Detector for 
Geostationary Earth Radiation Budget Applications,” M.S. Thesis, Virginia Polytechnic 
Institute and State University. 

Schmidt, W., 1995, “New Manufacturing Technology Improves Thermopile Sensors,” Laser 
Focus World, August. 

Sharma, A. and Cogley, A. C., 1982, “Radiative Heat Transfer in a Completely General 
Plane-Parallel Environment,” International Journal of Heat and Mass Transfer, Vol. 25, 
No. 4, pp. 523-534. 

Smith, S. M. and Howitt, R. V., 1986, “Survey of Material for an Infrared-Opaque Coating,” 
NASA Technical Memorandum 88204. 

Weckmann, S., 1997, “Dynamic Electrothermal Model of a Sputtered Thermopile Ther- 
mal Radiation Detector for Earth Radiation Budget Applications,” M.S. Thesis, Virginia 
Polytechnic Institute and State University. 



Vita 


Ira Joseph Sorensen was born^^^^B He attended Selma High 

School, from which he graduated in 1990. He received a Bachelor of Science in Mechanical 
Engineering from California Polytechnic State University, San Luis Obispo, where he was 
a member of Tau Beta Pi and an officer of ASME, in 1996. He relocated to Blacksburg, 
Virginia, in 1997 to begin his graduate studies at Virginia Polytechnic Institute and State 
University. Over the course of his graduate studies, he served as a research assistant in the 
Thermal Radiation Group, as well as a NASA Langley Researcher Summer Scholar (LARSS) 
and a NASA Graduate Student Researcher Program (GSRP) fellow. He received his Master 
of Science in Mechanical Engineering in December, 1998. He plans to continue his graduate 
studies in the Thermal Radiation Group to obtain his Ph.D. 


75 




